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Abstract: Cognitive Maps are a vital tool that can be used for knowledge representation and 
reasoning. Fuzzy Cognitive Maps (FCMs) are popular soft computing techniques used to model 
large and complex systems, and they can aid in explainable AI. FCMs, however, cannot model the 
indeterminacy that arises in a system due to various uncertainties. Neutrosophic Cognitive Maps 
(NCMs), upgraded FCMs that could model indeterminacy, was introduced to address this issue. 
NCMs are a generalization of FCMs, a field of cognitive science firmly based on neural networks. 
NCMs have been used to solve a wide range of problems. NCMs were introduced in 2002, and even 
after 20 years, NCMs do not have any supportive software or package, toolbox, or visualization 
software like FCMs. The main reason for the absence of dedicated software is due to the 
indeterminacy concept ‘I’ and how it has to be handled. This paper presents the dedicated Python 
package created for handling the functioning of NCMs. The modelling software presented in this 
paper aids in visualizing the NCMs as a signed digraph with indeterminacy that is directed signed 
neutrosophic graph. This package implements a sample case study using NCMs. 
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1. Introduction 


Fuzzy theory is a branch of mathematics that deals with vagueness and uncertainty in decision- 
making [1]. Fuzzy sets and logic model complex problems involving imprecise terms or partial truths. 
It has many applications in engineering fields, the healthcare sector, economics and social science, 
which pertain to real-world problems. Fuzzy logic and its models have many applications in various 
fields, such as engineering, artificial intelligence, medicine, economics, and social problems. It can 
help model multifaceted problems that involve human knowledge, preferences or emotions [2]. 

A Fuzzy Cognitive Map (FCM) [3] represents a mental landscape within which the connections 
between the nodes (e.g. events, concepts, resources or attributes) are used to compute the “strength 
of impact" of these elements. FCMs are signed fuzzy digraphs introduced by Bart Kosko [3]. 

FCMs have been used to analyze several socio-economic, healthcare and decision-making 
problems. The applications and extensions of FCMs are vast and widely researched; a few are 
presented here. 

In [4], the authors explore using FCMs as an agency for collective decision-making and how 
FCMs can capture the cognitive models and group beliefs of different stakeholders. FCMs were used 
as a learning assessment tool in [5] to understand the planning by children by stimulating cognitive 
function. 
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In [6] Multi-Agent Genetic Algorithm (MAGA) is proposed to optimize convergence error for 
learning FCMs. A GIS-dependent crisis management tool to predict earthquakes in Tehran using 
FCMs was proposed in [7]. 

An evolutionary algorithm, known as ]IBMTEA-FCM was proposed in [8] for learning large-scale 
FCMs. A qualitative analytical method using FCMs to specify causal-effect links of the 
interdependent SDGs by considering the long-term effect of COVID-19 was presented in [9]. 

In [10], the PRescriptiVe FCM (PRV-FCM) was introduced, based on FCMs and metaheuristic 
algorithms, to develop prescriptive models. In [11], three federated learning approaches were 
combined with FCMs for mortality prediction and treatment prescription with severe dengue cases. 

Iran’s population’s health was analyzed in [12] using FCMs. Since the concept of health is a 
complex and comprehensive system, other sector policies profoundly affect health. [13] proposes a 
decision support system for crowdfunding using FCMs. The various extensions of FCMs have been 
systemically reviewed by [14]. 

Several software packages are available for modelling FCMs. One such software, Mental 
Modeller , helps build FCMs intuitively and easily. After creating the models, decreasing or 
increasing the model’s elements allows us to examine various change tactics. FCMexpert is a software 
tool for FCM-based scenario analysis and pattern classification presented in [15]. Over 10 FCM 
extensions were handled by supporting interoperability in the FCM extensions in [16]. 

Python packages are also available for modelling FCMs. FCMpy [17] is a recently introduced 
open-source package for building and analyzing FCMs. Notably, FCMPy allows simulating the 
system behaviour using qualitative data to create fuzzy causal weights, applying ML algorithms to 
modify the FCMs matrix to aid in classification, and executing scenario examination by simulating 
theoretical interventions. 

The package also helps apply ML algorithms (e.g., nonlinear and active Hebbian learning, 
deterministic learning and genetic algorithms) to adjust the FCM weight matrix. 

Neutrosophy is a branch of philosophy investigating neutralities’ origin, nature, and scope and 
their interactions. Florentin Smarandache introduced neutrosophy in the 1990s [18]. Neutrosophy 
regards a proposition, hypothesis, concept, event, or entity depending on the modelling. 
Neutrosophy is the basis of the neutrosophic set, logic, probability, and statistics. Indeterminacy “I" 
is a concept in neutrosophy that measures the degree of neutrality or uncertainty of a proposition, 
event, theory, entity or concept. 

Neutrosophic Cognitive Maps (NCMs) are an extension of FCMs that can handle indeterminate 
relationships between two concepts, obtaining more significant sensitive results. It was introduced in 
[19] to analyse diverse social issues. NCMs have been modelled considerably on the AI focus to mimic 
the thinking-humanly approach. Here, it is unsupervised data and has a limited set of features. 

Over two decades, many investigators have utilized NCMs to analyze diverse problems like 
situation analysis [20], pest analysis [21], transgressions against people experiencing homelessness 
[22], and imaginative play in children [23]. FCMs and NCMs on COVID variants were compared in 
[24]. SWOT analysis and NCMs were combined to analyze organic farming in India [25]. [26] 
proposed triangular NCMs and uses them in multistage decision-making with a use case of 
evaluation. NCMs and cloud data were used in [27] to detect violence, and several datasets were 
used. NCMs and FCMs were compared in this analysis, and it clearly states that NCMs are better at 
handling indeterminacy than FCMs. 

Dynamic NCMs [28], enhanced cuckoo search and ensemble classifier were presented for 
acquiring the profile of the gene expression, differentiating between the individuals affected with 
rheumatoid arthritis and possible control subjects. [29] proposed a technique combining pest analysis 
based on fuzzy and neutrosophic logic to analyze the food industry. 
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In [30], the various factors impacting the paper-packaging industry are analyzed to provide a 
notional representation using NCMs since sustainable supply chains can be attained with repeated 
product use and recycling. 

In [31], NCMs were used to analyze the various causes and effects that lead to violent behaviour. 
NCMs were used to decide the elements that enable proper decision-making to confirm a precise 
diagnosis of conversion disorder [32]. The various factors that affect homeless people were analyzed 
using NCMs in [33]. A neutrosophic sociogram-based NCM approach was introduced in [34]. FCMs 
and NCMs have been used in health care to analyze dengue fever [35]. 

NCMs have been applied in various domains such as health, agriculture, engineering, social 
problems, business, law, environment, and medicine. The substantial advantage of NCMs over other 
cognitive maps is their capacity to capture data realistically and consistently represent expert 
opinions, making them a valuable tool for decision-making strategies where there is an advanced 
degree of indeterminacy or uncertainty. NCMs can be constructed either by the data-driven approach 
or experts’ opinion. 

NCMs have not been integrated with machine learning algorithms like FCMs. So, little research 
combines various machine learning algorithms that have been utilized in adjusting weights in an 
NCM, like in FCM. Despite the various applications of NCMs, there is no dedicated software or 
python package for NCMs. This paper presents a dedicated software and python package for 
constructing, analyzing and visualizing NCMs. 

The dedicated modules of our modelling software function in the following way: 

1. Generating the neutrosophic graph and the related connection matrix using expert opinion 

(a) from linguistic terms 
(b) from edge weights 
(c) from a file as input from the user. 

2. Visualizing the NCMs as a neutrosophic digraph. 

3. Simulating the NCMs using various state vectors. 

4, Analysis of various case scenarios for given NCMs. 


The paper is organized as follows: Section one is introductory. Section two recalls the working 
of NCMs and their construction. Section 3 provides the software-based visualization and simulation 
of NCM for a case study. The conclusions and results are presented in the last section, together with 
suggestions for future research. 


2. Working of Neutrosophic Cognitive Maps (NCMs) 


When data is unsupervised, and the association between two concepts is indeterminate, the 
indeterminacy can be captured by using neutrosophy [19] introduced the concept indeterminacy in 
FCMs, called NCMs. The basic properties of NCM are recalled to make this section self-contained. 

NCM is a digraph with concepts as nodes and their causal relationships as edges. These concepts 
can be events, strategies, or policies as nodes of the graph and relationships as edges, where each 
concept is mathematically represented as a neutrosophic vector from the neutrosophic vector space. 
Every node, in its vector form, is represented by (x1, , Xn); x; € {0,1,/}, where 0 is off state, 1 is on 
state and I is the indeterminate state. 

Consider two nodes Ng and N, of the NCM; their relationship is given by the edge eg,. Every 
weighted edge e,, is from {-1, 0, 1, J, where 0 means no impact, positive value means increase in 
Nq implies increase in N,, similarly decrease in N, implies decreases in Ny. If eg, takes a negative 
value like —1, it implies that an increase in Ng implies a decrease in N,, or similarly, a decrease in 
Na implies an increase in N,. The edge weight is assigned a value / if the effect from N, to Np, can 
not be determined. The edge weights of simple NCMs are from {-1, 0, 1, I}. The NCM’s adjacency 
matrix is denoted by N(E) = (e ap), Cay € {-1,0,1,]} . 
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NCMs are expert opinions constructed based on data obtained from the expert, where they 
identify the factors/concepts relevant to the domain and associated causal relationships in terms of 
numbers and indeterminacy or linguistic terms. 

Generally, the neutrosophic connection matrix is obtained directly from the digraph of the 
expert’s opinion. Here, we have introduced the concept of dealing with neutrosophic linguistic terms 
using the following logic. Algorithm 1 provides neutrosophic edge weights from the linguistic terms. 


Algorithm 1: Generate NCM Matrix 


Data: matrix «+ Empty Matrix; row <- Empty List; n < No of nodes 
Result: NCM Matrix; 
i,j — 0 while i# ndo 
string < “” 
while 7 # n do 
term < input linguistic term 
if i= j then 
| element + 0 
else 
if term = ”-VH” then 
| element 
else if term 
| element 
else if term 
| element 
else if term 
| element 


+ random(—5, —2.75) 
= "-H” then 
+ random(—2.75, —2) 
= "-M” then 
< random(—2, —1.5) 
= ”-L” then 
+ random(—1.5, —1) 
else if term = ”-VL” then 
| element < random(—1,0) 
else if term = ”"NC” then 
| element <0 
else if term = "+VH” then 
| element < random(2.75, 5) 
else if term = "+H” then 
| element + random(2, 2.75) 
else if term = "+M” then 
| element < random(1.5, 2) 
- 
- 


else if term "+L" then 

| element random(1, 1.5) 
else if term "+VL” then 

| element random(0, 1) 
else if term = "NaN” then 


| element <— I 
string < string + eleme nt+',’ 
j¢jti 
for « € temp do 
ife2=”"I" then 
| continue 


row <— row.append(temp) 
i¢it+1 
iO 
while i #4 n do 
matrix < matrix.row_insert(i, row{i]) 
'—#+1 
return matrix 


In the case of linguistic terms, these edges can be considered as negative very high (-VH), 
negative high (-H), negative medium (-M), negative low (-L), and negative very low (-VL). It implies 
a negative influence when the decrease of the influence of the node N, results in the increase of the 
influence of the node N, or the increase of N, results in the decrease of N,. Similarly, these edges 
can be considered as positive very high (+VH), positive high (+H), positive medium (+M), positive 
low (+L), and positive very low (+VL). It implies a positive influence when the decrease of N, results 
in the decrease of N, or the increase of N, results in the increase of Nj. If the linguistic term is no 
causality, 0 is used; if it is indeterminate, “NaN" (not a number) is used. These linguistic terms are 
converted to numerical values using simple assignments as follows: 

e If term = “-VH" then the element gets a random value from [-5, -2.75). 
e If term = “-H", then the element gets a random value from [-2.75, -2). 
e If term = “-M", then the element gets a random value from [-2, -1.5). 
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elf term = “-L", then the element gets a random value from [-1.5, -1). 
elf term = “-VL", then the element gets a random value from [-1, 0). 

e If term = “-NC" then the element gets 0. 

e If term = “+VH" then the element gets a random value from (2.75,5]. 
elf term = “+H", then the element gets a random value from (2, 2.75]. 
elf term = “+M", then the element gets a random value from (1.5,2]. 
elf term = “+L", then the element gets a random value from (1, 1.5]. 
elf term = “+VL", then the element gets a random value from (0, 1]. 

e If term = “NaN" then the element gets I. 

The values generated by the algorithm range from [-5, 5]; since the non-indeterminate edge 
weights of NCM are from [-1, 1], the edge weights are normalized before the NCM/neutrosophic 
digraph is constructed. 

The expert opinion obtained can be used to generate the neutrosophic adjacency matrix of the 
NCMs. The generation can be done by using linguistic terms. According to the linguistic terms, the 
neutrosophic matrix can be obtained using random values between the ranges. 


Example 1 Consider the graph of the NCM, using the seven concepts (or nodes or attributes) 
Ly, L2,..,L7 and the expert opinion is obtained in terms of linguistic terms, using Algorithm 1, the 
weights for the edges are provided. 

For illustration, assuming that the connection from L, to L7 has a very high positive influence, 
and the connection between from L, to L3 is indeterminate, and the connection from L, to L, isa 
very high negative influence. The edge weights generated by the algorithm using the linguistic terms 
are tabulated in the Table 1: 


Table 1: The edge weights assigned by the algorithm 


Ly 
Lo 
Ly 
Ly 
Ls 
Le 
Ly 


The NCM will be given by Figure 1. Each edge is weighted and directed. The dashed lines are used 
to represent indeterminate edges. The weight of each edge is given in the graph. 


Figure 1: An illustration of the neutrosophic directed graph 
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It is then normalized to weigh between [—1,1] and indeterminacy I. The connection matrix for 
the neutrosophic directed graph is given in Equation 1. 


0 -05 J 03 O 0.7 O 
0 0 0 O 0.3 O 1 
0 0 0 O 1 0.4 0 
0 I 0 0 0 0 —0.9 
Enormalized = 0 0 0 0 0 0 -02 (1) 
0.1 O 0 O 0 0 I 
0 0 0 O 0 0 0 


The edge weights can also be obtained from the expert. The edges of simple NCMs are from {-1, 0, 
1, B. 


Example 2 Consider a simple NCM given by an expert with seven concepts Ly, L2,...,L7 as the nodes of the 
directed neutrosophic graph. The expert opinion is obtained in terms of edge weights. 


Figure 2: Neutrosophic directed graph for simple NCM 


The edge weights are from {—1,0,1,/}, and the indeterminate edges are represented by dotted lines. 
The NCM’s connection matrix is denoted by N(E) =(e gy), where eg, € {0,1, —1, I}. 


0 -1 r 1 0 YT O 
0 0 00 1 0 1 
0 1. 1h OA Ae 
oF 0 0 0 0 =1 

ae 0 0 0 0 0 0 -1 (2) 
11 0 0 0 0 71 
0 0 1 0 0 0 0 


The notion of state vector, dynamical system and its functioning are described. The neutrosophic 
state vector S = (sj,...,5,);5; € {0,1, 1}; where 0 indicates the off state, 1 is in the on state, and I 
indicates the indeterminate state. Let LL, LjL3, L3Ll4, .. , Lalp be the NCM’s directed edges. 
Given the edges creating a cycle, the NCM becomes cyclic; otherwise, it is acyclic. Consider L,L3, 
Tolgusuplaln, tobe cyclic, if node L,_ is on the influence will flow via the existing cycle and cause 
L, to be on again. This state of equilibrium of the dynamical system is called as a hidden pattern. 
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Various state vectors with different nodes in on state are considered to activate the system 
Consider the NCM with feedback given in Figure 2; its neutrosophic adjacency matrix is N(E) given 
in Eq. (2). 

The state vector S, = (1, 0, 0,, 0) where L, is in on state is considered. The data needs to be 
transformed by N(E), so we multiply 5, by N(E). 

The state vector multiplied by the neutrosophic matrix N(E) is given in Algorithm 3. In the 
given algorithm, S is a 1-dimensional matrix (row vector) and the neutrosophic adjacency matrix 
N(E) is a 2-dimensional matrix denoted by B. The resultant vector Res of the multiplication of 
S x N(E) is returned. 

An example illustrates this: Consider the graph in Figure 2 and its related connection matrix. 
Take the state vector S; = (1000000). The state vector S,; is multiplied with the neutrosophic 
adjacency matrix N(E). 


0 <1 7 1 8 1 6 

0 0 00 10 1 

00 0 0 FT 4 0 

= 0 © 0 0 0 0 -1 
S,xN(E)=(1000000)x), 9 9 9 9 9 <1 (3) 

110 0 0 0 7 

0 0 1 0 0 0 0 
=(0 —111010) (4) 


Algorithm 2: Multiply 
Data: 1D Matrix S and B < 2D Neutrosophic Adjacency Matrix N(E) 
Result: Resultant Matrix 
Res < 1D Matrix 
r < S.size 
R+« B.size 
i,j <0 
while i#r do 
while j 4 R do 
Res{i] — Res{i] + (S{i] * Blj][{i)) 
j¢jtl 
iti+l 
return Res 
After obtaining the resultant vector Res, it must undergo the threshold and update operations. By the 
definition of NCM, the threshold and update operation is denoted by the © symbol. 
It is important to note here that working with J (indeterminate) needs to be done carefully. As 


by the definition of J 
Ixl=P=!I (5) 


Any power of I gets mapped to J, as shown in Algorithm 3. 


Algorithm 3: UpdateIPower 


Data: V < state vector 
Result: V state vector after updating IJ’s 
4<0 
forieV & jE€bdo 
| ifi=J? then ic lI 
return V 


Similarly, during the updating and threshold operation, any constant into J is also mapped into 
I, as shown in Eq. (6) 
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nxl=1 (6) 


The resultant vector from the multiplication of the state vector with N(E) is thresholded and 
updated. 
Let X= S ,N(E) =(s 1,5 2,,8 ) is thresholded by replacing s; accordingly to the equation. 


1 if Si >t 
s,=40 if s;<t (tis a suitable positive integer) (7) 
I if s;not an integer 


If the resultant vector is 1+ I, it is made to 1, whereas in case | — 1, it is thresholded and updated to 
0. The resultant X is updated to ensure that the state considered on in the initial state vector S, is on 
in the resulting vector. Here, it is updated to make the concept N , as 1 in the resulting vector. 


The algorithm for thresholding and updating is given in Algorithm 4. 


Algorithm 4: ThresholdandUpdate 
Data: X < resultant vector and ¢ threshold value, Initial On state 
Result: Thresholded and updated resultant vector 
X + updateIPower(X) 
len + X.size 
iO 
while i 4 len do 
temp-expr < X [i] 
if temp_ecpr = J then 
| temp.expr + 0 
if temp_ezpr > t then 
| X[ij<t 
else if temp_expr = 0 then 
if X{i] 40 then 
| Xfi} J 
else 
| X[ij)] <0 
t<-it+l1 
X (state — 1) + 1 
return X 


Here, the Algorithm 4 illustrates where only one node is taken in the on state in the initial state 
vector. 

The thresholding and updating operation is mathematically denoted by ©. For example, 
consider the resultant vector of S, x N(E), the thresholding and updating result in S,. 


5, X N(E) = (0 —111010) 6 (10/1010) =S, (8) 


Here, during the threshold operation, —1 is made 0, and during the update operation, the very 
first state is made on again. 

Generally, in any NCM, more than one concept/ node can be considered in the on state. This will 
deal with the combined effect of both states being on. The proposed model can handle more than one 
node in the on state, and it can work to analyse the effect of a combination of various nodes. The 
multiplication of resultant vectors with the neutrosophic adjacency matrix N(F) will continue until 
the resultant vector yields a fixed point or limit cycle. 

If the NCM settles to a neutrosophic state vector repeating in the form 


Si GS S2 GS Sj ig Si S Si, 


Kandasamy et al., NCMPy: A Modelling Software for Neutrosophic Cognitive Maps based on Python Package 


Neutrosophic Systems with Applications, Vol. 13, 2024 9 
An International Journal on Informatics, Decision Science, Intelligent Systems Applications 


Then the dynamic system’s equilibrium is called NCM’s limit cycle. Suppose, S,N(E) © S2 
(where © denotes the resultant vector of $;N(E) whichis thresholded and updated) and for S ,N(E) 
we repeat the same procedure until we attain the fixed point / limit cycle. 

Two vectors are compared in Algorithm 5. 


Algorithm 5: Compare 


Data: a,b - 1D vectors 
Result: Boolean value 
ij <0 
forica & j€bdo 
if i~j then 
| return false 
return true 


A fixed point or limit cycle is attained when both the vectors under comparison are the same. 
For example, if we compare 5S, with S,, we can see that S, # S3. 


(1000000) #(10/1010)=S, (9) 
Since the fixed point nor limit cycle is reached, i.e., S, # Sz, the process is continued. 
S,x N(E) = (11111141 -14+D°(1111110) =S3 (10) 


Since the fixed point still needs to be reached, i.e., S3 # Sz, the process is continued. 


S,xX N(E)=(1+/ 20 21.1 2214+7 -1+21D 
6(1111110)=S, (11) 


Here S3 = S,. The fixed point has been reached. 


The Algorithm 6 is used to determine if a limit cycle is reached. It compares with the previous 
resultant vectors using the previously described in compare Algorithm 5 to check if the limit cycle is 
reached. In case it is reached, it returns a true or a false. 


Algorithm 6: Check_cycle 
Data: b < list of vectors 
Result: Boolean value 
n < b.size 
i<O 
while i #4 n—1do 

jeitl 

while j 4 n do 
vl < Di] 
v2 < d{j] 
if compare(vl, v2) = true then 

| return true 

jejtl 

itit+l 

return false 


To illustrate an example of limit cycle, consider the state vector P, = (0100000), and the 
connection matrix N(E) as given in Eq. (2), that is the related adjacency matrix for the NCM given in 
Figure 2. 
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P, x N(E) = (0000101) 6 (0100101) =P, 
P, x N(E) = (0010100) 6 (0110100) =P, 
P; X N(E) = (0000210) 6(0100110)=P, 
P, x N(E) =(1100101 9(110010ND =P, 
P; x N(E) = (0-12*11110) 
6 (0111110) =P, 
P, x N(E) =(11+100/4+11I-1) 
6 (1100110) =P, 
P, x N(E) = (11-111111?) 
6 (1111111 =P, 
Pg X N(E) = (1122 +11 1+12+*I1-1) 
6 (1111110) =P, 
Py»X N(E)= (1 +1-111141141 2-1 
6 (1111110) =P, 
Py X N(E) = (112 11+12*11-1) 
6 (1111110) =P, = Py (12) 


Since P,, = Py, the iteration is stopped since a limit cycle has been achieved, enabling the 
determination of the hidden pattern. The limit cycle is as follows: Pg gives Pg, Py gives Pio, Pio 
gives P,,; that is same as Pp. 

The process of multiplication of the resultant vector with the matrix N(E£) is continued until a 
limit cycle / fixed point is reached. 

The Algorithm 7 takes an adjacency/connection matrix and state vector (an integer as input 
denoting the state to be activated) and an integer denoting the threshold values as parameters. 


Algorithm 7: Iteration 
Data: an Adjacency matrix FE, state < state to be activated, t threshold value 
Result: list of vectors 
len + E. size 
cl ¢ list{len] filled with 0 
cl[state — 1] + 1 
start + Matriz(cl) 
flag — false 
vectors + Empty List 
while flag = false do 
y + multiply(start, F) 
y <- thresholdAndUpdate(y, t) 
vectors < vectors.append(y) 
start + y 
flag < check_cycle(vectors) 
return vectors 


The Algorithm 7 continues multiplying the resultant vector with the N(E) until a fixed point 
or a limit cycle is reached. It is dependent on the several previously described algorithms. Every 
concept must be made in the active state to capture the hidden pattern and understand the effect of 
the concept/ node on others to analyse the NCM thoroughly. Only in NCMs can we signify that the 
clout of a node on different nodes can be indeterminate, and this vision needs to be revised in the 
case of FCMs. 

For the given example, the results have been tabulated in Table 2 for various state vectors. 
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Table 2: Determining the Hidden Pattern 


Input state vector limit cycle / fixed point Resultant vector 
(1000000) fixed point (1111110) 
(0100000) limit cycle (1111110) 

(1111110) 
(1111110) 
(1111110) 

(0010000) fixed point (1111110) 
(0001000) fixed point (01011100) 
(0000100) fixed point (0000100) 
(0000010) fixed point (1111110) 
(0000001) fixed point (1111111) 


It is seen that only the active state of concept Lz results in a limit cycle—the rest results in a fixed 
point. 

The conclusions that can be drawn from the equilibrium state of the dynamic system are as 
follows: When node L, is active, all nodes are either indeterminate or on state. When L, is in the on 
state, it results in a limited cycle, which affects all nodes other than L7. Similarly, when nodes 
L3,L4,L¢ or L7 alone are in the on state, other nodes are either indeterminate or on state. Whereas 
when node L. is in on state, no other node is affected; all of them remain in the off state. 


3. Description of the modelling package 


The overall flow of the modelling software is as given in Figure 3. The first module is for the 
input module, which can either be linguistic term-based or edge-weight-based. 


Causality in Linguistic 
Terms 
« Assign Edge weights 
according to terms 


Creating the 


Input No of Neutrosophic A scr Simulate the 
Nodes Graph using the saath dynamic system Results 
edge Weights 


State vector 


i e \Svs ened ve + S2=S1*N(E) 
, > 1000000 * Threshold and update 


¢ Direct Edge Weights a phe 6 08 operation 
from Expert opinion oo71r0 ¢ Limit Cycle/Fixed point 


Figure 3. Various modules of NCM modelling software. 


1. Input module: There are two methods in which the model can be created using expert opinion. 
They can enter by either method as described below: 

- Linguistic terms: The user can enter the linguistic terms to describe the relationship between 
two concepts in NCM as very negative or positive. According to the previously discussed 
algorithm, the edge weights are assigned based on the linguistic term. 

- Edge weights: The user can instead directly enter the edge weight to denote the causality 
between two nodes. 

2. Neutrosophic digraph: The neutrosophic digraph is generated and visualised using the edge 
weights obtained from the user. 
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3. Connection matrix: The related connection matrix N(E) is obtained from the neutrosophic 
digraph of the NCM. 

4. Dynamical system: Using various state vectors, the dynamical system is simulated to analyse 
the effect of the on state of the various nodes. 

5. Results: The effect of the on state of various nodes and combination of various states is 
consolidated. 

The NCMpy python package performs a simulation of the dynamic system—a detailed 
description of the package is given in the next subsection. 


3.1 Description of the NCMPy Python Package 
The flowchart of the NCMpy package is given in the following figure. 


eapue State Vector updateIPower () CSTON 
Fixed Point/ 
A 4 A Resultant Resultant 
iteration () ——> multiply () — li —» thresholdAndUpdate () —— ae —» check_cycle() —yYes—> Limit Cycle 


| 7 | 


Figure 4: Flowchart of the NCM package 


The SymPy library has been used to handle indeterminate values in modelling the NCMs in this 
proposed package. In SymPy package, the names E, I, O, 5, N, and Q collide with names already 
defined in the package. Hence, J can not be used to represent indeterminacy, so J is used instead of 
I. Throughout the coding snippets, J is used instead of I. 

The various functions/modules are described here. 


1. compare(x, y): This function takes two vectors as input and checks for their equality; it is 
dependent on Algorithm 6. The code snippet is as follows: 
def compare(x, y): 
res = True 
for k in range(len(x)): 
af xiky) t= yifikiy: 
return False 
return True 


Assume that compare function takes S$; = (111110) and Sg = (III 1110) as x and y. In 
that case, compare(x, y) will return True, in case it was S; = (1/11/11) and Ss =CUIIIII0) 
under consideration, then it would return False. 


2. check_cycle(b): This function is dependent on Algorithm 6; it is used to check if the resultant 
vector is a fixed point or limit cycle. 
def check_cycle(b) : 
for i in range(len(b) - 1) : 
for j in range(i + 1 , len(b)) : 


v1 = b[i] 
v2 = b[j] 
if(compare(v1 , v2) == True) : 


return True 
return False 
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The fixed point is achieved when the resultant vector is the same as the previous resultant vector, 
that is P;N(E) © P;. The limit cycle is achieved when the recently calculated state vector is the same 
as any one of the previously calculated resultant vector, which results in a cycle. 


P:N(E) & Pisa; PitaN(E) © Pita) PN (E) & Pi 


In Example 2, considering the active state P; = (0100000) results ina limit cycle. 


3. updateIPower(x): This function takes a 1D vector as input and converts the indeterminate 
quadratic polynomial to a linear polynomial. This function is based on Algorithm 3. This function is 
used by the threshold and update function. 
def updateIPower(x): 

for i in range((np.shape(x))[1]): 
if x[i].find(3**2): 
x[i] = x[i].xreplace({J**2: J}) 
return x 


Consider the resultant vector; 
R, = CU *#*2JJ +110 ** 2 2/) 
For a sample scenario. The updateIPower(R_1) will change this vector P, into 


R, = (JJ +110) 2) 


4. thresholdAndUpdate(X,threshold_value,state) : This function takes a 1D vector and a 
threshold value as a parameter and updates each vector value according to the defined thresholding 
operation. Also, the updation operation must see to it that the node which was on in the initial state 
is on in the next state and so on in the resultant state also; if not, it is set to 1 again. 


def thresholdAndUpdate(X , threshold_value, state) : 
X=updatelIPower (X) 
for i in range((np.shape(X)[1])): 
if (X[i].find(-1)): 
x[{i]=0e 
temp_expr=X[i].subs(J,@) 
if (temp_expr>=threshold_value): 
X[i]=X[i].subs(X[i],threshold_value) 
elif (temp_expr==0): 
if(X[i].find(J)): 
X[{i]=3 
else: 
X[{i]=X[i].subs(X[i],9) 
if isinstance(state, list): 


activated_states = [i for i, x in enumerate(state) if x == 1] 
for s in activated_states: 
X[s] =1 
else: 
X[state - 1] =1 
return X 


Xp XN(E)= (I+ 17 +1111 412*«11-1)9 (1111110) =X, 
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5. iteration(E,state,threshold_value) : This function takes an adjacency/connection matrix, an 
integer denoting the state to be activated and an integer denoting the threshold values as parameters. 
This function is based on the Algorithm 7 . 


def iteration(E , state, threshold_value = 1) : 
if isinstance(state, list): 
start = sym.Matrix(state) 
else: 
cl = np.zeros((np.shape(E)[1])) 
cl[state - 1] =1 
start = sym.Matrix(c1) 


flag = False 

start = start.T 

vectors = [] 

while flag == False : 
y = multiply(start , E) 
# Performing the thresholding operation on output vector y 
y = thresholdAndUpdate(y , threshold_value, state) 
vectors.append(y) 
# Updating start vector to start with new state vector 
start = y 
# Checking for cycle among state vectors 
flag = check_cycle(vectors) 

return vectors 


Consider Example 2, where the working out is done with active state P, = (0100000); it 
iterates until the limit cycle is achieved. 
3.2 Modelling NCMs for sample case study 


The opening page for the modelling software is given in Figure 5. Here, the user can select the 
option of working with linguistic terms or directly entering the neutrosophic edge weights given by 
the expert as shown in Figure 6. 


WE HELP YOU TO 


Model Your 
NCMs 


Figure 5: Homepage 


WHATWEDO —- ABOUT US 


How would you like to give input? 


Version 112 


Figure 6: Option Selection 
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3.3 NCMs using Linguistic Terms 


The neutrosophic linguistic terms are obtained from the user, as shown in Figure 7. 


Inguisti 


npu c terms 
ee Add Nodes 


Figure 7: Sample input page for the linguistic terms related to the neutrosophic directed graph 


The NCMpy package runs through the Alogrithm GenerateNCM 2 and creates the necessary 
edge weights and normalises them, providing the neutrosophic bigraphs. The neutrosophic-directed 
graph of the NCM is given in Figure 13. 


Table 3: The edge weights assigned 


Li is L; iy ie be io 
Li 0 4.88 I 3.58 1.91 -2.39  -2.81 
i 4.22 0 I 2.13 I 2.74  -4.68 
i. I I 0 -4.77 I 3.23. 1.93 
i” 2.98 1.58 -2.41 0 241 | 181 | 3.11 
ie 123 I I 42 0 2.12 I 


Le -1.20 -4.36 4.03 1.71 -2.40 0 2.85 
Ly -2.03  -4.77 1.34 3.89 I 4.44 0 


The resultant neutrosophic connection matrix of the graph is 


[ 

[@, 0.98, J, 0.72, 0.38, -0.48, -0.56], 
[0.84, 0, J, 0.43, J, -0.55, -.94], 

[J, J, @, -0.95, J, 0.65, 0.39], 

[@.60, 0.32, -0.48, 0, 0.48, 0.36, 0.62], 
[0.24, J, J, 0.84, @, -0.42, J], 

[-0.24, -0.87, 0.81, 0.34, -0.48, @, 0.57], 
[-0.41, -0.95, 0.27, 0.78, J, 0.89, Q] 
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0.39 


0.24 


Figure 8: Neutrosophic directed graph 


To show the plotting capacity of the visualizing module, we have taken a matrix with all connections 
for the sample. The threshold value is obtained from the user. According to the threshold value set 
by the user, the thresholding and updating of state vectors are done. 


Figure 9: Insert the threshold value and state vector 


The threshold value of 1 was taken here, and the following results were obtained for various 
state vectors. The working out for the first state vector S$, = (1000000), where the concept C, is 
on is shown here: 
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FOR ACTIVE STATE 1 

Matrix([[@, 0.98, 1.0*J, @.72, 0.38, -@.48, -@.56]]) 

Matrix([[0, @.98, 1.0*J, @.72, 0.38, -@.48, -0.56]]) 

Matrix([[@, 0.98, 1.0*J, 0.72, 0.38, -@.48, -@.56]]) 

Matrix([[@, @, 1.0*J, 0.72, 0.38, -@.48, -@.56]]) 

Matrix([[@, ®, J, 0.72, @.38, -@.48, -@.56]]) 

Matrix([[@, @, J, @, 0.38, -@.48, -0.56]]) 

Matrix([[0, ®, J, @, @, -@.48, -0.56]]) 

Matrix([[@, ®, J, @, @, @, -@.56]]) 

After Thresholding and Updating 

Matrix([[1, ®, J, @, @, @, @]]) 

Matrix([[J**2, J**2 + @.98, J, 0.72 - @.95*J, J**2 + 0.38, 0.65*J - 0.48, @.39*J - @.56]]) 
Matrix([[J, J + 0.98, J, @.72 - @.95*J, J + 0.38, 0.65*J - 0.48, @.39*J - @.56]]) 
Matrix([[J, J + 0.98, J, @.72 - @.95*J, J + 0.38, @.65*J - 0.48, @.39*] - 0.56]]) 
Matrix([[J, ®, J, @.72 - @.95*J, J + 0.38, @.65*J - 0.48, 0.39*J - 0.56]]) 
Matrix([[J, @, J, 0.72 - @.95*J, J + 0.38, @.65*J - 0.48, 0.39*J - @.56]]) 
Matrix([[J, @, J, 0, J] + 0.38, @.65*J - 0.48, 0.39*J - 0.56]]) 

Matrix([[J, ®, J, @, @, @.65*J - 0.48, 0.39*J - @.56]]) 

Matrix([[J, 0, J, @, @, @, @.39*J - 0.56]]) 

After Thresholding and Updating 

Matrix([[1, ®, J, @, @, @, @]]) 

[Matrix([[1, @, J, @, @, @, 0]]), Matrix([[1, 0, J, @, @, 9, @]])] 


Similarly, the working out for each and every state vector is carried out. For state vector B, = 
(0 100000), where the concept Cz is on, the resultant vectors will be: 


FOR ACTIVE STATE 2 

Matrix([[0.84, 0, 1.0*J, 0.43, 1.0*J, -0.55, -@.94]]) 

After Thresholding and Updating 

Matrix([[@, 1, J, @, J, ®, @]]) 

Matrix([[J**2 + @.24*] + 0.84, 2*J**2, J**2 + J, 0.43 - @.11*J, J**2 + J, 0.23*J - 0.55, J**2 + @.39*J - 0.94]]) 
After Thresholding and Updating 

Matrix([[@, 1, J, @, J, @, @]]) 

Resultant vector(s) [Matrix([[@, 1, J, ®, J, 0, @]]), Matrix([[@, 1, J, 0, J, @, @]])] 


For state vector G,; = (0010000), where the concept C3 is on, the resultant vectors will be: 


FOR ACTIVE STATE 3 

Matrix([[1.0*J, 1.0*J, @, -0.95, 1.0*J, 0.65, @.39]]) 

After Thresholding and Updating 

Matrix([[J, J, 1, @, J, @, @]]) 

Matrix([[2.3*9, J**2 + 2.0*J,. 3*3**2,. 2.0*3 = 9:95, J**2 + 1,473, 0.65 = 1,473, JF*2 = 1.5*3 + 0.39] ]) 
After Thresholding and Updating 

Matrix([[J, J, 1, @ J, @, @]]) 

Resultant vector(s) [Matrix([[J, J, 1, @, J, 0, @]]), Matrix([[J, J, 1, ®, J, @, 0]])] 


For state vector X,; = (0001000), where the concept C, is on, the resultant vectors will be: 


FOR ACTIVE STATE 4 

Matrix([[@.60, 0.32, -@.48, @, 0.48, 0.36, 0.62]]) 

After Thresholding and Updating 

Matrix([[@, @, @, 1, @, @, Q]]) 

Matrix([[@.60, @.32, -0.48, @, 0.48, 0.36, 0.62]]) 

After Thresholding and Updating 

Matrix([[@, @, @, 1, @, @, @]]) 

Resultant vector(s) [Matrix([[@, 9, @, 1, 0, @, 0]]), Matrix([[@, 0, 0, 1, 0, 0, @]])] 


For state vector Y, = (0000100), where the concept C; is on, the resultant vectors will be: 


FOR ACTIVE STATE 5 

Matrix([[@.24, 1.0*J, 1.0*J, 0.84, 0, -@.42, 1.0*J]]) 

After Thresholding and Updating 

Matrix([[@, J, J, ®, 1, ®, J]]) 

Matrix([[J**2 + 0.43*] + 0.24, J**2 + @.05*J, J**2 + 1.3¥*J, 0.26*J + 0.84, 3*J**2, @.99*J - 0.42, @.45*J]]) 
After Thresholding and Updating 

Matrix([[@, J, J, @, 1, @, J]]) 

Resultant vector(s) [Matrix([[®, J, J, @, 1, @, J]]), Matrix([[®, J, J, @, 1, ®, J]])] 


For state vector Z, = (0000010), where the concept C, is on, the resultant vectors will be: 


FOR ACTIVE STATE 6 

Matrix([[-@.24, -0.87, 0.81, 0.34, -0.48, @, 0.57]]) 

After Thresholding and Updating 

Matrix([[@, ®, @, @, @, 1, @]]) 

Matrix([[-@.24, -0.87, 0.81, 0.34, -0.48, @, 0.57]]) 

After Thresholding and Updating 

Matrix([[@, @, @, @, @, 1, @]]) 

Resultant vector(s) [Matrix([[@, ®, @, 0, @, 1, 0]]), Matrix([[®, 0, @, @, 0, 1, @]])] 
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FOR ACTIVE STATE 7 
Matrix([[-0.41, -@.95, 0.27, 0.78, 1.0*J, 0.89, @]]) 
After Thresholding and Updating 
Matrix([[@, @, @, @, J, @, 1]]) 
Matrix([[0.24*] - @.41, J**2 - 0.95, J**2 + 0.27, 0.84*J + 0.78, J, 0.89 - 0.42*J, J**2]]) 
After Thresholding and Updating 
Matrix([[@, @, ®, @, J, ®, 1]]) 
Resultant vector(s) [Matrix([[@, 0, ®, @, J, @, 1]]), Matrix([[@, 0, ®, @, J, @, 1]])] 
For state vector A; = (000000 1), where the concept C; is on, the resultant vectors will be: 
The result vector in each case has been shown. Results regarding the resultant vectors that can 
be discussed 
1. Maximum Influence: The nodes C2, C3 and Cs; are the most influential since they do 
affect many other nodes and turn to an indeterminate state. 
2. Least Influential nodes: The nodes C, and Cg are the least influential since they do not 


affect any other node than itself. 
3.4 NCMs using Edge Weights 
The edge weights are obtained from the expert, as shown in Figure 10. 


Input Neutrosophic Matrix 


J0j0000 


1900000 


“1500500 


1000000 
0000s0-1 


-1900000 


Figure 10: Input screen 


Using the edge weights obtained, the NCM is created. It is the visualization of the same represented 
as a NCMs digraph as shown in Figure 11. 


Figure 11: Visualization of the neutrosophic graph 


The related neutrosophic connection matrix N(E) is given in Equation 2. The obtained matrix is 


[@, -1, 'J", 1, @, 1, @], 
[@, 9, @, @, 1, @, 1], 
(0, ®, ®, @, 1, 1, @], 
[@, ‘J’, 0, 0, @, @, -1], 
(0, ®, @, @, 0, @, -1], 
(1, 1, @, @, @, 9, ‘J'], 
(9, @, 1, @, 0, @, @] 
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The threshold value of 1 was taken here, and the following results were obtained for various 
state vectors. The working out for the first state vector S$, = (1000000), where the concept C, is 


on is shown here: 


FOR ACTIVE STATE 1 


Matrix([[@, -1.e@@e@e@eeeeeee, 1.0*3J, 1.e@@@@eeGeGeG00, 8, 1.eQeeee0eE00000, 0]]) 


After Thresholding and Updating Matrix([[1, ®, J, 1, @, 1, @]]) 
Matrix([{[1, J, J, 1, 3, 3 +1, J - 1])) 
After Thresholding and Updating Matrix([[1, J, J, 1, J, 1, @]]) 
Matrix([[1, 3, J, 1, 2*3, J +1, J - 1))) 
After Thresholding and Updating Matrix([[1, J, J, 1, J, 1, @]]) 


Resultant vector(s) [Matrix ([[4,, 0, 3, 1,: O,, 2, @]))}; Matrix(€[[4,; J, 3, 2, Jy 15. O9]), Matrix([[t,. 3, Jy 2, 35 25 OJ1)) 
The resultant vectors for all active states are given in the snippet below: 

FOR ACTIVE STATE 2 

Resultant vector(s) [Matrix([[@, 1, 0, @, 1, 0, 1]]), Matrix([[@, 1, 1, ©, 1, 0, @]]), Matrix([[e, 1, e, 

®, 1, 1, @]]), Matrix([[1, 1, @, @, 1, @, J]]), Matrix([[@, 1, J, 1, 1, 1, @]]), Matrix([[1, 1, 0, ®, 1, 

J, @]]), Matrix([[J, 1, J, 1, 1, 1, 3]]), Matrix([[1, 1, J, J, 1, J, @]]), Matrix([[J, 1, J, 1, 1, 1, 

@]]), Matrix([[1, 1, J, J, 1, J, @]])] 

FOR ACTIVE STATE 3 

Resultant vector(s) [Matrix([[@, @, 1, @, 1, 1, @]]), Matrix([[1, 1, 1, @, 1, 1, @]]), Matrix([[1, @, 1, 

1, 1, 1, J]]), Matrix([[1, J, 1, 1, 1, 1, @]]), Matrix([[1, J, 1, 1, 1, 1, 6]])] 

FOR ACTIVE STATE 4 

Resultant vector(s) [Matrix([[@, J, @, 1, 0, @, @]]), Matrix([[@, J, ©, 1, J, @, @]]), Matrix([[e, J, 2, 

1, J, @, 6]))] 

FOR ACTIVE STATE 5 

Resultant vector(s) [Matrix([[@, ®, 0, ®, 1, 0, @]]), Matrix([[@, ®, @, ®, 1, @, @]])] 

FOR ACTIVE STATE 6 

Resultant vector(s) [Matrix([[1, 1, @, @, 9, 1, J]]), Matrix([[1, @, J, 1, 1, 1, 1]]), Matrix([[1, J, 1, 

41, J, 4; OT); Matrix([[t, 3, 35. Dy By. 1, 09]),, Matrix([(d,. J. J, 2, 3; 2, 0]]), Matrix€[ fa, 3, 3; 2, Js 

1, @]])] 

FOR ACTIVE STATE 7 

Resultant vector(s) [Matrix([[@, @, 1, @, @, @, 1]]), Matrix([[@, @, 1, @, 1, 1, 1]]), Matrix([[1, 1, 1, 

®, 1, 1, 1]]), Matrix([[1, ®, 1, 1, 1, 1, 1]]), Matrix([[1, J, 1, 1, 1, 1, 1]]), Matrix([[1, J, 1, 1, 1, 


1, 1]))] 


The combination of two concepts being in the on state is also handled by the package done there. 
Similarly, several other combinations can be done, according to the user’s choice. A sample of the 
combinations of various concepts begins in on state, and the working is given in the snippet. 


Activating state: [1, 1, 0, 8, @, 0, 9] 


(Matrix([[1, 1, J, 1, 1, 1, 1)]), Matrix([[1, 1, 1, 1, 1, 1, @]]), Matrix([[1, 1, J, 1, 1, 1, @]]), 
Matrix([[4, 1, 3, 1, 1,1, 0]])] 

Activating state: [1, 9, 1, 0, @, @, 9] 

[Matrix([[1, @, 1, 1, 1, 1, @]]), Matrix([[1, J, 1, 1, 1, 1, @]]), Matrix([[1, J, 1, 1, 1, 1, @]])] 
Activating state: [1, 0, 0, 1, @, 0, @] 

[Matrix([[1, @, J, 1, @, 1, @]]), Matrix([[1, J, J, 1, J, 1, @]]), Matrix([[1, J, J, 1, J, 1, @]])] 
Activating state: [1, 90, 0, @, 1, @, 9] 

[Matrix([[1, ©, J, 1, 1, 1, @]]), Matrix([[1, J, J, 1, 1, 1, @]]), Matrix([[1, J, J, 1, 1, 1, @]])J 
Activating state: [1, 0, 0, 0, @, 1, 9] 

[Matrix([[1, ®, J, 1, ®, 1, J]]), Matrix([[1, J, J, 1, J, 1, @]]), Matrix([[1, J, J, 1, J, 1, @]])] 
Activating state: [1, 0, 0, 0, @, @, 1] 

[Matrix([[1, @, 1, 1, @, 1, 1]]), Matrix([[1, J, 1, 1, 1, 1, 1]]), Matrix([[1, J, 1, 1, 1, 1, 1]])] 
Activating state: [@, 1, 1, 0, @, @, 9] 

[Matrix([[@, 1, 1, @, 1, 1, 1]]), Matrix([[1, 1, 1, ©, 1, 1, J]]), Matrix([[1, 1, 1, 1, 1, 1, 3]]), 


Matrix([[1, 1, 1, 1, 1, 1, @]]), Matrix([[1, 1, 1, 1, 1, 1, @]])] 


The complete NCMPy python package will be made available online. 


4. Conclusions and discussions 


As the world is moving towards no code or less coding paradigms, a python package for NCM 
will aid and help mathematicians, social scientists, economists, strategists and other policymakers in 
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analyzing real-world problems without worrying about the mathematical background or 
computational complexities of NCM. A visualization tool and python package to help in the working 
of NCM were presented in this paper. The NCMpy package and modelling software provide the 
functions essential to study problems involving indeterminacy, which can be done using NCMs. 

This package and modelling tool is open-source, written in python, straightforward to 
implement, and provides the required functionality for handling models with indeterminacy. 

This tool implementation is a collaboration with the founding and leading experts from the field 
of NCMs. This tool will facilitate research and enable new researchers and scientists to apply NCMs 
in their projects which involve indeterminacy. We plan to update our library and constantly welcome 
all scientific community contributions. 
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